Electronic structure of the a-phase in paramagnetic Fe-V alloys. Experimental and 

theoretical study. 
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The electronic structure of cr-phase Feioo-xVi compounds with 33.3 < a; < 60.0 was calculated 
from the charge self-consistent Korringa-Kohn-Rostoker method. For the first time, charge densities 
Pa(0) and electric field gradients were determined at Fe nuclei, that occupy five nonequivalent lattice 
sites. The highest pa(0) values were found on sites A and D, and the lowest one on site B, the 
difference ranging between 0.f62 and 0.f74 s-like electrons per Fe atom for x — 33.3 and x — 60, 
respectively. The calculated quantities combined with experimentally determined site occupancies 
were successfully applied to analyze ^^Fe Mossbauer spectra recorded on a series of 8 samples in a 
paramagnetic state. 



Studies of complex intermetallic phases, where several 
nonequivalent crystallographic sites exist, are of a great 
interest Q&Sli- Especially challen ging in this context 
are disordered alloy systems with the complex structure 
such as cr-phases. The latter, whose archetypal exam- 
ple is that in the Fe-Cr system belongs to an im- 
portant class of tetragonal close-packed crystallographic 
structures Q viz. the Frank-Kasper phases 0. They 
are characterized by a high coordination number (CN), 
varying from CN=12 to 15 in the cr-phase. The unit cell 
(space group PA2/mnm) contains 30 atoms distributed 
in a non-stoichiometric way over 5 crystallographic sites, 
commonly labeled A, B, C, D and E. The importance 
of this class of phases is further reinforced by the fact, 
that they exhibit topological properties similar to sim- 
ple metallic glasses due to an icosahedral local arrange- 
ment [1]. Consequently, they can also be regarded as 
very good approximants for dodecagonal quasicrystals 
Q. Indeed, close similarity between vibrational prop- 
erties of a one-component cr-phase on one hand and a 
one-component glass with the icosahedron local arrange- 
ment on the other, was recently reported [3]. 

The complex structure and a lack of stoichiometry 
also cause the interpretation of experimental results ob- 
tained for cr-phase samples to be a difficult task. Here, 
a cr-phase magnetism which can be found in some alloy 
systems like Fe-Cr and Fe-V is of a particular interest 
[3, [lo| . In order to understand it properly, one should 
know sublattice magnetic properties. An access to the 
latter is experimentally possible with microscopic tech- 
niques such as Mossbauer Spectroscopy (MS) or Nuclear 
Magnetic Resonance (NMR) . We are not aware of using 
NMR in studying Fe-containing cr-phases, but MS has 
been already applied and distinction between electronic 
properties of nonequivalent sites was not obtained due to 
low resolution in the measured spectra d, However, 
backed by dedicated electronic structure calculations, rel- 
evant information on Fe charge-densities and electric field 
gradients (EFG) on all nonequivalent sites has been al- 
ready provided in the case of the cr-Fe-Cr compounds in 



a paramagnetic state 

In this Letter we report results obtained with similar 
calculations for the cr-phase in Feioo-^Va: alloys, also in 
the paramagnetic state. This series of compounds is es- 
pecially well-suited for verifying electronic structure com- 
putations, since a composition range of the Fe-V cr-phase 
existence is several times wider than that in the Fe-Cr 
system. The non-spin polarized KKR results have been 
used to analyze ^"^Fe-site Mossbauer spectra recorded at 
295 K (paramagnetic state) in the cr-phase compounds 
Feioo~a;V;j; with 34.4 < a; < 59.0. 

The Mossbauer spectrum associated with each site can 
be characterized by the following spectral parameters: 
amplitude /, line width G, isomer shift (center of grav- 
ity), IS, that measures a charge-density at nucleus of the 
probe atom, and quadrupole splitting, QS, that gives in- 
formation on the electric field gradient (EFG) in the lat- 
tice. Taking into account one parameter for a background 
and the fact that in the unit cell of the sigma there are 
five sites (sublattices) one needs 21 parameters to de- 
scribe the overall spectrum. However, as the ^^Fe-site 
Mossbauer spectrum of the cr-phase in the paramagnetic 
state has no well- resolved structure [lo| , unique determi- 
nation of these parameters from the spectrum analysis 
itself appears impossible and additional information is 
required. The latter can be partly obtained e. g. from 
a neutron diffraction experiment that has not only con- 
firmed the presence of the probe ^^Fe atoms on all 5 
sites but also yielded precise knowledge on their rela- 
tive population on these sites I4I . Assuming the Lamb- 
Mossbauer factor to be site independent, the relative con- 
tribution of each subspectrum ascribed to the site should 
be equal to the corresponding values determined from the 
neutron experiment. Setting the latter as constraints, re- 
duces the number of free parameters from 21 to 17 which 
is still to high to allow the unique refinement of the spec- 
trum. Theoretical calculations, as those presented else- 
where, are of a great help in this respect |llj . 

It seems that the KKR method combined with the 
coherent potential approximation (CPA) would be best 
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FIG. 1: (Online color) Fe-site charge-density, pa(0), for five crystallographic sites versus the number of NN-Fe atoms for (a) 
Feee.rVsa.a, (b) Fe4oV6o and (c) Fe54Cr46. The lines stand for the best linear fits to the data that represent pyi(0)-values for 
particular atomic NN-configurations. The average value of pa(0) for each site is indicated by a horizontal bar on the right-hand 
axis on which are marked corresponding /5- values in mm/s. 



adapted to investigate electronic structure of chemically 
disordered cr-phase as it occurs in the Fe-V system. 
However, such self-consistent calculations are very time- 
consuming due to the presence of Fe/V atoms disorder 
on all sites as well as a high multiplicity (8i) of three 
sites giving rise to a long process of CPA self-consistency 
cycles. In addition, such results may only give us the 
quantities averaged over all possible atomic configura- 
tions. This can be important when studying theoret- 
ically magnetic properties of these compounds, but it 
appears rather useless in the determination of the spec- 
tral parameters in the paramagnetic state. We have pro- 
posed an alternative approach to face the problem [ll|, 
in which the electronic structure of the cr-phase was cal- 
culated in a reasonably large number of ordered approx- 
imants, namely atomic configurations. These configura- 
tions should represent the most probable Fe and V ar- 
rangements in the disordered cr-phase unit cell and cover 



wide range of the Fe-nearest neighbour (NN) numbers 
but account for the fact that some configurations are 
much less probable than others [ill] . Consequently, we 
have taken into accout those NN configurations whose 
overall probability was greater or equal 0.95. Such prob- 
ability cut-off has permitted to significantly reduce the 
configurations number and diminish computation time 
to a reasonable level. In practice, all calculations were 
carried out using the lowest symmetry simple tetrago- 
nal unit cell (PI), in which the atoms had the same po- 
sitions as in the original symmetry group (P4:2/mnm), 
but they were occupied either by Fe or V atoms. In 
the case of the Feioo-xV^: alloys, where the cr-phase ex- 
ists for 33 < a; < 60, two border compounds have been 
chosen, i.e. Fe2oVio corresponding to x = 33.3, and 
Fei2Vi8 equivalent to x = 60. In order to obtain the 
most probable atomic configurations that satisfy the con- 
dition on the nearest neighborhood appearing with the 
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overall probability > 0.95, it was enough to take into ac- 
count 26 and 17 different configurations for the x = 33.3 
and a; = 60 compounds, respectively. More details on the 
calculations can be found elsewhere 

Fig. 1 presents the charge-density at ^^Fe nuclei, 
/9a(0), and the corresponding isomer shift, IS, calculated 
for the specified cr-phase compositions for nonequivalent 
sites. Analyzing all dependences one can come to the 
following conclusions: (i) pa{^) is characteristic of each 
site and it linearly decreases with the Fe-NN number at 
a rate typical for each site and composition, (ii) the high- 
est charge-density have the sites A and D, and the lowest 
one the site B, and the difference between them depends 
on the composition but hardly on the number of Fe-NN 
atoms, (iii) the value of the charge-density at a given site 
depends on the alloy composition. The decrease of the 
charge-density with the increase of the Fe-NN atoms ob- 
served for all sites and compositions means that there is 
an effective charge transfer from V atoms to Fe atoms. A 
similar effect was observed for disordered bcc Fe-V alloys 
[l^ . Using a fitted scaling constant from 14 1, one can 
correlate the calculated charge-densities with the num- 
ber of s-like electrons. Thus the average charge density 
at Fe nuclei on the site A is equal to 0.175 for x — 33.3 
and to 0.205 s-like lectrons for a; = 60. For site B the 
corresponding values are 0.013 and 0.031. It is also in- 
teresting to estimate in this way the difference in the 
charge-density between the extreme sites i.e. A and B. 
By doing so one arrives at 0.162 for x = 33.3 and 0.174 
for X = 60 i.e. the difference hardly depends on the com- 
position. 

The charge-densities for any other x-value lying be- 
tween 33.3 and 60 can be determined as a weighted aver- 
age of the parameters obtained for the two border com- 
positions. Taking into account /^-values calculated for 
different compositions, atom configurations and sites one 
can calculate an average IS- value, < IS >i {i — A, B, 
C, D, E) for each nonequivalent site. The < IS >i values 
obtained in this way are marked on the right-hand scale 
in Fig. 1 as horizontal bars. The probability of each site. 
Pi as well as the aforementioned < IS >i allow deter- 
mining the average isomer shift for selected composition 
< IS >, 



<IS >=^P,< IS >i 



(1) 



As illustrated in Fig. 2 theoretically calculated < IS >- 
values are in a good agreement with the measured center 
shift, < CS >. The data presented here also give a clear 
evidence that substituting Fe by V in the cr-phase struc- 
ture has quite similar effect on the Fe-site charge- deiisity 
values as already observed in the a-phase structure [13| . 
where it also increased the Fe- site charge-density. 

The good agreement between the experimental and 
calculated data give us a confidence that our electronic 
structure calculations deliver reliable description of the 
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FIG. 2: (Online color) Average center shift, < CS > versus 
vanadium concentration, x, as measured (circles) and as cal- 
culated (solid line). The dashed line indicates the best-fit to 
the measured values. 



hyperfine parameters in the cr-phase samples. One may 
expect that a more detailed information obtained for the 
sublattice charge-density is also credible. 
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FIG. 3: (Online color) Quadrupole splitting, QS, as deter- 
mined for each site and vanadium concentration, x, from the 
analysis of the measured spectra with the protocol described 
in the text. Values obtained previously for a cr-Fe-Cr com- 
pound are indicated as reversely-filled symbols for compari- 
son. 

Concerning the second spectral parameter i.e. QS, it 
was calculated based on an extended point charge model 
as described in detail elsewhere jTlj. The obtained QS- 
values for each site and composition are shown in Fig. 3. 
It is clear from the figure that also QS is characteristic of 
each site, maybe except for the sites B and D for which 
QS"- values are very close to each other whatever the con- 
centration X. Another interesting feature is that except 
for the A site the QS'-values remain constant vs. com- 
position for all sites. The greatest QS-value was found 
for site E ('^0.45 mm/s) and the lowest one for C (~0.20 
mm/s). For comparison, QS'- values obtained previously 
for the cr-Fe-Cr compound 11[ are marked by reversely- 
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filled symbols in Fig. 3. Interestingly, for the sites A, C 
and E they match very well the values calculated in the 
CT-Fe-V, while for the sites B and D they are significantly 
smaller and differ one from another. Finally, using for 
the relative intensities of the five subspectra the values 
as found from the neutron diffraction experiment [l^ on 
one hand, and taking for IS and QS the values obtained 
from the present calculations on the other hand, we have 
successfully refined the measured ^''Fe Mossbauer spec- 
tra. In the fitting procedure, each subspectrum was re- 
garded as composed of double lines with the same QS 
and /^-values linearly dependent on the number of the 
Fe-NN atoms and probabilities obtained from the bino- 
mial distribution. It should be mentioned that only five 
free parameters were needed for the analysis, i.e. back- 
ground, total spectral area, IS for site B (to adjust the 
refined spectrum to the used source of the gamma rays), 
line width and a proportionality factor between QS and 
an energy shift due to the quadrupole interactions 
The agreement between the experimental and calculated 
spectra is excellent for all measured compositions of the 
cr-Fe-V system. 
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FIG. 4: (Online color) ^^Fe-site Mossbauer spectra recorded 
at 295 K on two samples with (a) x = 37 and (b) x = 46. The 
best-fit spectrum and five subspectra are indicated by solid 
lines. 

In summary, our results of the KKR electronic struc- 
ture calculations obtained for the cr-Fe-V system lead to 
conclusion that the proposed procedure, based on the hy- 
perfine parameters calculated from the first principles in 



combination with the experimental data for the atom dis- 
tribution over the five sublattices, gives information on 
the charge-density and the electric field gradient in the 
complex disordered alloy system. It should be mentioned 
here, that this information was obtained for all five sub- 
lattices separatelly, taking into account the real atomic 
concentrations on them. 

The progress and adventage of our calculations in com- 
parison to those carried out by others, in which homoge- 
neously occupied sublattices were considered, is that we 
have taken into account much more various atomic con- 
figurations that has enabled us to get more deep, precise 
and realistic insight into the electronic structure of the 
CT-phase. Our calculations have been, in addition, verified 
by a comparison with experimental data. A very good 
agreement achieved can be taken as evidence that the 
procedure we have applied gives a good approach to the 
real electronic structure of the studied system. It can be 
also of a practical importance while analyzing ill-resolved 
Mossbauer spectra recorded on samples with a complex 
structure. 
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